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NTIS  ~~ GRAJfel 
DTIC  TAB 
Unannounced  □ 

Justification _ 


By— - 

Distribution/ 

Availability  Codes 
(Avail  and/or 
Dist  j  Special 


notations : 


TL  stands  for  the  set  of  all  integers, 
K  for  the  real  numbers, 

1?0  :=  {n  €  TL  ,  n  >  0), 

:=  {n  e  TL  ,  n  >  l), 
x Ay  :=  min(x,y}, 

xV  y  :=  max(x,y} , 

(x)+  :=  max(0,x}, 

(x)  :=  min(0,x}, 

fx"]  :=  m±n{i  e  TL  |i  >  x}, 

[xj  :=  maxti  €  TL  |i  <  x). 


for  all  ne»i;  (*)  :=  1;  (*)  :=  0  for  all  z  {  WQ, 


For  the  values  of  a  function  V:  •>  IR  we  use  both  notations  v(i) 

and 


Tables  of  Significance  Points  for  the  Variance-Weighted 
Kolmogorov-Smirnov  Statistics 
By 

Heinrich  Niederhausen 


1 .  Introduct ion . 


Let  xi»***»^  be  i.i.d.  random  variables  with  continuous  distri¬ 
bution  function  F  and  empirical  distribution  function 

M 


F  (x)  -  M"1  l  1 


i-1 


,  ,(X.)  . 

(-<*>, x]  i 


The  goodness-of -f it  statistic 


Fx(x)-F(x) 


__  *  aw 

W  «  sup  - - 


has  been  shown  to  be  asymptotically  minimax  (with  respect  to  a  certain 

loss  function)  by  A. A.  Borokov  and  N.M.  Sycheva  (1968).  They  also  give 

some  exact  significance  points  and  the  asymptotic  distribution  of 

/M  Ww.  Beside  Ww,  we  consider  the  following  related  statistics: 

M  M 


\f  (x)-F(x)| 

WM  -  sup  . . ----- 

01<F(x)<92  /f(x)  (1-F(xj  ) 


Fy(x) -F(x) 

WM  -  sup  ,  ' 

9L<F  (x)<02  /F(x)  (l-F(x)) 
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If  (x)-f(x)! 

WM  "  SUP  - -  .  ,  - 

61<Fx(x)<92  »^F(x)  (l-F(x)) 

+  |Fx(x)-Fy(x) | 

J  =  sup  . -  -  ■  ZZZ 

’  e1£Fv(x)<e2  /fv(x)(i-fv(x» 

|fx(x)-Fy(x) | 

N  =  sup  ~7=H=ZZZ=Z 

011Fv(x)<02  /Fv(x)(1-Fv(x>) 


Fx(x)-Fy(x) 

Slip  - - 

91-FX(x)-62  ^v(x)(l-Fv(x)) 


|Fx(x)-Fy(x) 

J  =  sup  .... -  , 

*  01<Fx(x)<02  /fv(x)(1-Fv(x)) 


where  Y^,...,YN  is  a  second  independent  sample  with  the  same  distribu¬ 
tion  function,  and  V  ,  is  the  combined  sample.  We  call  all 

these  statistics  variance-weighted  Kolmogorov-Smirnov  tests.  In  [10], 
we  derived  some  methods  to  compute  the  exact  distribution  of  such  tests. 
Using  those  methods,  we  computed  tables  for  the  significance  points  of 


wj,  JH  wM,  w+,  vfi  wM, 


WM,N’ 


WM,N  and 


Let  Z  be  any  of  the  eight  statistics  in  (1).  Let 


P(z)  -  P(Z  <  z)  . 


For  each  a  *  .9,  .95  and  .99  we  try  to  find  z  such  that  P(z  )  ■  a. 

a  a 

But  the  variance-weighted  Kolmogorov-Smirnov  distributions  are  discontinuous, 
even  in  the  one-sample  case.  Therefore,  we  give  P(z  )  and  z  in  the 

\X  \X 

tables,  where  P(z  )  is  smaller  than  a.  After  each  z  ,  a  single 
“Cl 

digit  D  is  printed.  If  the  last  digit  of  z^  is  increased  by  D,  a 
z  is  obtained,  such  that  P(z  )  >  a.  P(z  )  is  also  listed.  All  numbers 
are  rounded  in  the  last  digit. 

In  all  the  tables  we  chose  0^  =  ®  for  9  “  0,  0.01,  0.05, 

4*  4* 

0.1  and  0.25.  In  Ww,  W„  and  Ww  ,,  we  have  to  take  the  supremum 

M  M  M,N 

over  0  £  Fx(x)  £  1-0.  Thus,  we  have  to  replace  0  by  d/M,  where  the 
integer  d  is  chosen  such  that  d/M  comes  close  to  0  (see  (3.2)). 
Analogously,  replace  0  by  d/(M4-N)  in  w"*”  and  W  . 

M  y  N  n  y  N 

In  the  two  sample  case,  all  tables  are  given  for 
M  *  2,3,4,. ..,10  ;  N  -  2,3,4, ...,M 

M  =  15, 20, 25,..., 50;  N  «  M,M-l,M-2, . . . ,M-5 

M  =  100,500  ;  N  =  M  . 

If  for  small  M  the  table  for  a  certain  0  does  not  differ  from  the 
preceding  table  (with  smaller  0),  then  this  part  of  the  table  is  omitted. 

In  the  one  sample,  the  eame  values  of  M  are  used,  but  the  sample 
length  M  *  500  is  omitted.  The  computer  proved  to  be  too  slow  for  this 
case  (and  the  desired  accuracy) . 
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of  a  two  sided 


For  large  sample  sizes,  a  significance  value  z^ 
statistic  can  be  approximated  by  z^  ^  the  corresponding  one  sided 
statistic.  The  larger  the  9,  the  better  the  approximation.  Despite 
"bad"  asymptotic  behavior,  this  approximation  is  practically  satisfying 
even  for  9=0.  For  this  case,  the  computer  drawing  on  the  next  page 
shows  z  g  and  z+g$  for  M  =  2 . 80,  where 

M"5  <  z.9>e„0  *  -9  “d  p<^  <  <  <95)e-o  *  -95  • 

We  demonstrate  on  page  6  what  happens  if  z+g$  (which  is  much  faster  to 

compute)  is  used  as  an  approximation  for  z  Q:  We  plotted  M  against 

P(*^i  W  <  z+ne'tn-n  for  M  “  2,..., 100.  To  illustrate  the  effect  of  a 

. _  00 

larger  9,  we  plotted  also  P(vM  W„  <  z  )  .  The  jumps  come  from 

m  —  .yo  0=1 

the  approximation  of  9  by  d/M. 

We  repeat  that  part  of  [10],  which  is  necessary  to  understand  the 
algorithms.  Chapter  3  gives  an  overview  over  the  significance  points  by 
small  tables.  The  large  tables  are  computed  in  the  same  way. 

All  computations  are  done  on  a  pdp  11  computer,  using  16  signifi¬ 
cant  digits.  The  reader  can  compare  the  tables  for  ✓S  W*  with  table  1 
and  2  of  A. A.  Borokov  and  N.M.  Sycheva  (1968).  With  the  expection  of  one 
printing  error,  their  numbers  differ  at  most  by  1  in  the  last  given  digit. 
The  case  8=0  for  W  has  been  considered  by  M.  Noe  (1972).  His 
method  of  computation  is  close  to  ours  for  all  two-sided  one  sample 
statistics  in  (1),  except  for  special  cases,  see  3^. 
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for  all  n  6  M,  .  Obviously,  f  (v  )  =  6A  .  hence,  (f  )  has  roots 
x  n  n  Vjii  n 

in  v  (see  (A. 14))- 

Denote  by  the  1-th  order  statistic  of  a  size  M  random  sample 

from  U(0,1).  If  u  .  <  1,  then 

M-l 

(1.1)  ^(Mm-i^  =  p^vi  -  U(i)  -  Mi-1  Vi  = 

=  ’ 
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1.2.  Recursions. 


For  this  section  we  assume  y^  -  1  without  loss  of  generality. 

Ti"“  k 

With  *  x  / (n-k) !  in  (A. 13),  the  algorithm  A.l  was  found 

by  M.  Noe  (1972).  Given  that  y^  is  constant  for  k  =  0,...,K,  say, 

it  may  be  possible  to  use  an  explicit  formula  for  f  (y  )  (i=0,...,K). 

i  K. 

The  same  is  true  for  the  values  p^  in  the  following  application  of 
(A. 16):  Define  p^  :=  1  and 

a.?)  p^-T  (‘)<-i)1-k-1<Pk-v.)‘-k  pk  . 

Then 


P(v.  £  U  £  for  all  i  =  1, .  ...M)  =  pM  . 


Alternatively,  we  get  from  corollary  A. 2: 


pM  -  detCc/.Hy.  1-v.)f‘J+1).  .  M 

M  j-1  j-1  i  +  i,j=l,...,M 


V.A.  Epanechnikov  (1968)  found  recursion  (1.2)  and  G.P.  Steck  (1971) 
independently  derived  this  determinant  and  many  applications.  See  also 
E.J.G.  Pitman  (1972)  for  another  proof. 


Remark:  Depending  on  the  accuracy  of  the  computer,  (1.2)  should  be  used 
only  for  small  M  because  of  the  alternating  summation.  In  algorithm  A.l 
the  summation  does  not  alternate,  but  compared  with  (1.2)  the  amount  of 


computations  is  approximately  squared!  In  the  computation  of  signifi¬ 
cance  points  it  often  occurs  that  £  p^.  The  same  is  true,  of  course, 

in  both  one  sided  cases.  Theorem  A. 2  (with  i  :■  0)  yields  for  any  real 
function  a  on 


(1.3)  Pj  "  J!tj>o(0j)  "  JQ  (k)(aj-Uk)J"kpk  for  a11  jml . M  * 

From  £  p^  we  see  that  t^  ^  equals  for  j  »  0, . . .  ,M  the  Sheffer 
sequence  for  D  with  roots  in  v.  Given  that  a  p^  ^  for  all 
j  ■  1,...,M,  the  summation  in  (1.3)  is  non-alternating.  But  how  to 
compute  tj  q(Oj)?  In  the  simple  case  v  =  0  (one  sided  tests)  we  get 
j!t^  Q^j)  =  a j  -  With  0^  :=  Pj  ^  Steck’s  formula  (1971,  (2.3))  is 
obtained.  See  the  previous  section  for  other  closed  forms.  We  suggest 
the  following  procedure  for  general  v  (with  £  Pq)  r 

Choose  0=1.  Thus,  for  all  j  =  1,...,M, 

j!tjyi)  -  P(vi  <  U(i)  for  all  i  -  1, . . . ,  j)  =  j !f jj) (1)  , 

If  (f^)  is  the  p^-Sheffer  sequence  for  D  with  roots  in  0,  where 
n 

■  l-v(j-i)  for  all  i  «■  0, . . . , j  .  Hence,  (1.3)  can  be  applied  to 
compute  t  _(1)  (we  choose  again  0  =  1): 

3>" 

p£3)  -  1  and  p<j)  -  1  -  T  (^(j-k)1"^5  for  all  i«l,...,j  . 

k=0 

Thus,  j!tj  0(1)  -  pj33  for  all  j  ■  1,...,M.  Finally,  enter  again 
(1.3)  and  compute  pM  from  Pq  ”  1  and 

P1  -  pfJ)  ~  l  (k)(1_Wk)J"kpk  f°r  a11  3  “  1 . M  ’ 

J  J  k“0 
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1.3.  Renyl  type  distributions. 

In  applications,  the  test  distributions  seldom  occur  in  the  form 
of  (1.1).  But  if  our  method  is  applicable  at  all,  they  are  easily  trans¬ 
formed  so  that  one  of  the  following  two  lemmas  can  be  used: 

Lemma  1.1:  Let  f  and  g  be  monotone  non-decreasing  functions  from 
[0,1)  into  itself  such  that  f  <  g  and 

f(0)  <  a/M  <  b/M  <  g(l)  **  1 

for  two  fixed  integers  a  and  b.  Then 

P(f (Fy(x))  <  x  <  gCFyCx))  Va/M  <  Fu(x)  <  b/M) 


=  P(vi 

.  i  V)  ± 

Wi-1  V1 

-  1 . M) 

1 

r  0 

for  all 

i  »  0, . . . ,a-l 

1 

•H 

.  f(i/M) 

for  all 

i  *  a, . . . ,b 

1 

[  f (b/M) 

for  all 

i  >  b. 

g(a/M) 

for  all 

i  -  0, . . . , a— 1 

g(i/M) 

for  all 

i  »  a, . . . ,b 

1 

for  all 

1  >  b  . 

The  proof  is  obvious.  The  situation  in  the  following  lemma  is  much  more 
complicated. 
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Lemma  1.2:  Replace  a/M  and  b/M  in  lemma  1.1  by  any  real  ct  and  0 
such  that  0  £  a  <  0  £  1.  Then,  under  the  same  assumptions  about 

f  and  g , 


PCfCFyCx))  £  x  £  gCFyCx))  Va  £  x  £  0) 

-  P(v.  <  \J.  .  <  y  for  all  i  =  1,, 
x  —  (i)  —  i-1 


.,M) 


if 


0  for  all  i  *  0, . . .  ,ot 

f(i/M)  for  all  i  =  a^+l,...,6^ 

0  for  all  1  >  B, 


and 


for  all  i  *  0,...,a^-l 
g(i/K)  for  all  i  =  a  ,...,8  -1 


■  gU, 
k  1 


for  all  i  >  0 
“  g 


where 


af  “  max{k  £  M|f (k/M)  £  a) ,  0,  *»  max[k  £  M|f  (k/M)  £  0} 


a  =*  minCk  £  0|g(k/M)  £  a)  ,  0  =  min(k  £  o|g(k/M)  £  0}  . 

8  8 


(M'i 

Proof.  Denote  by  [0,1]  the  set  of  all  monotone  non-decreasing 

M 

ordered  vectors  u^ [0,1]  : 


u  =  (Uj,, . . .  ,uH>  such  that  0£u1£**'£uM£l  . 

(M) 


Define  the  subset  A  of  [0,1] 


by 


A  :*  {f (i/M)  £  x  holds  for  all  i  -  and  xttu^u^pftfa.B]} 

(uQ  :=*  0,  :=  1).  Then 

A  ■  (f(i/M)  £x  holds  for  all  1  **  a^+l,...,M  and  x€[u^,ui+^)fl[a,|3]} 

=  (f(i/M)  £  holds  for  all  i  =  a^+l,...,M  such  that  £  B) 

■  >  B,  and  f(i/M)  £  holds  for  all  i  =  ot^+1, ..  .,3^  such  that  £  B) 

“  {ug  +1  >  B,  and  f(i/M)  £  u^  holds  for  all  i * ocf+l, . . . ,B^}  . 

By  interchanging  the  roles  of  f  and  g  it  follows  analogously  that 


B  :=  {x  £  g(i/M)  holds  for  all  i=0,...,M  and  xe[ui>ui+1)f|[a,B]} 


(u^  <  a,  and  ui  £  g((i-l)/M)  for  all  i=  a  +1,...,B  }. 


P(AfiB)  ■  P(f(Fy(x))  £  x  £  g(Fy(x))  Va  £  x  £  B)  finishes  the  proof. 


Remark:  If  <  for  all  i**  0,...,M-1  in  the  lemmas  above,  look 
for  the  beat  applicable  method  in  1,2.  The  probability  is  zero 
otherwise. 
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2. 


Two  sample  tests. 


2.1.  Sheffer  polynomials  for  V. 

Denote  by  the  set  of  all  vectors  T  consisting  of  exactly 

i  ones  and  j  zeros.  For  each  T  =  (T^ . . -.T^)  e  define  the 

path  T^  of  T  by  :  =  0  and  T^  :=  Tk  for  all  J l  =  l,...,i+j. 

The  set  Gfd.j)  is  closely  related  to  empirical  distribution  functions: 

Let  Y^,...,Yn  be  M+N  continuous  and  i.i.d.  random  variables. 

Denote  the  monotone  non-decreasing  ordered  combined  sample  by  V-,...,V 

X  M+N 

Define  a.e. 


(2.1) 


if 

V*“Xi 

for  some 

1, 

1 

li£ 

M 

if 

VA  =  Yj 

for  some 

j. 

1 

<  i  1 

N  . 

Then  Tj  =  MF^CV^)  and  £-1^  =  NF^(V^).  Let  y  and  v  be  integer 
valued  function  on  155^  ,  -1  =  £  y^  and 


(2.2)  Vi-1-1  -Vi<  '‘i  l-^i  for  ®  • 

Then  fi(j)  =  #3'(i,j|v(T£)  <  *-T£  <  y(T£)  for  all  *-  0,...,i+j),  if 
(fn)  is  the  y-Sheffer  sequence  (with  variables  in  20  for  the  back¬ 
wards  difference  operator  V  (see  (A. 6))  with  roots  in  V  .  Hence, 

M+N  ~ ^ 

(2.3)  P(v(Tj)  <  fc-Tj  <  y(Tp  for  all  l  -  0 . M+N)  -  (  M  )  f^(N)  . 


13 


2.2.  Recursions. 


We  assume  y(M)  **  N  throughout  this  section.  From  the  definition 
of  a  y-Sheffer  sequence  (f^)  for  7  with  roots  in  V  we  get  the  following 
two-dimensional  recursion 


(2» 


fi<j)  -  | 


f  i(j-l)+f  ^(j) 
0 


for  all 


else  , 


\<  ilv, 


with  initial  values  fg(j)  =  1  for  all  j  £  y^,  and  f  (v^)  *  0  for 
all  i  e  11^  .  On  a  computer  with  unlimited  integer  precision,  this 
algorithm  may  be  slow  but  absolutely  accurate! 

The  one-dimensional  recursion  (A. 16)  is  left  to  the  reader.  From 
corollary  A. 2  one  gets  the  determinantal  solution 


P(v(t£)  <  £-T^  £  y(Tp  for  all  l  =  0,...,M4fl) 

L 


M+N 
V  N  1 


-1 


det 


'(Y 

i-J+l 


This  determinant  has  been  found  independently  by  G.  Rreweras  (1965)  and 
G.P.  Steck  (1969).  See  also  S.G.  Mohanty  (1971)  and  E.J.G.  Pitman  (1972) 
for  other  proofs. 

A  close  look  on  v  and  y  may  save  some  recursion  steps.  If 
v(M)  <  y(0)  the  outside  method  allows  non-alternating  summation  as 
described  in  1.2. 
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2.5.  Renyi  type  distributions. 

Lemma  2.1.  Define  f,  g,  a  and  b  as  in  lemma  1.1.  Then 


P(f(Fx(x))  <  Fv(x)  <  g(Fx(x))for  all  a/M  <  Fx(x)  <  b/H) 


v 


i 


P(v(Tp  <  l-T'z  <  p(Tp 
if 

f  _1 

4  f"(M+-N)f(i/M)*]-i-l 


for  all  2,=  0,...,M+-N), 

for  all  i  ■  0,...,a-l 
for  all  i  *■  a, . . . ,b 
for  all  i  >  b  , 


and 


f  “a 

for 

all 

i  » 

0 . a-1 

«  S  l_(MfN)g(i/M)j-i 

for 

all 

i  « 

a  y . . .  f  b 

1  N 

for 

all 

i  > 

b  . 

The  proof  is  obvious  from  2 . 1  . 

Lemma  2.2.  Define  f,  g,  a  and  b  as  in  lemma  1.1,  and  a^,  8^»  a 

3  as  in  lemma  1.2  with  o  :=  a/(M+-N)  and  3  :■  b/(>HN).  Then 
8 


P(f(Fx(x))  <  Fy(x)  <  g(Fx(x)>  for  all  a/(K+N)  <  Fy(x)  <  b/(M+N)) 

*=  P(v(Tp  <  l  <  y(Tp  for  all  2  «=  0 . M4tf)  , 

if 


-1 


for  all  i  ■  0, , 


.,a# 


|lMfN)f(i/M)"]-i-l 
1  b-0f-l 


for  all  i  *  a^+1, . . .  ,3  ^ 
for  all  i  >  3r  , 


and 
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and 


y 


i 


a-a 

g 

-  [_(H+N)s(i/M)J-i 

N 


for  all  i  •  0, _ ,ot  -1 

g 

for  all  i  -  a  .... ,3  -1 
g  g 

for  all  i  >  3 
~  g 


The  proof  follows  the  same  pattern  as  the  proof  of  lemma  1.2  and  Is  there¬ 
fore  omitted. 

Remark:  It  may  happen  that  v  or  y  in  lemma  2.1  or  2.2  violates  the 
monotonicity  conditions  (2.2).  In  this  case  define  the  "monotone  hulls" 

*•  A 

v  and  y  by 


(2.5) 


°o  :=  -1 

V.  :=  max{v.,v.  ,}  for  all  i  =  1.....M  , 
i  i  i-1 


and 


(2.6) 


yM  :=  N 

y^  :=  min{y^,yi+^}  for  all  i  «•  0,...,M-1  . 


A 

If  <  for  i  "  0,...,M-1,  look  for  the  best  applicable 

method  in  2 ■ 2.  The  probability  is  zero  otherwise. 
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5.  The  variance-weighted.  Kolmogorv-Smirnov  tests. 


We  defined  W.,  in  the  introduction.  Let 
M 


+  2i+s+[s2+4si(l-i)]1/2 

h  (i)  =  - 


2(l+s) 


and 


c  (y)  -  M(y  +[sy(1-Y)3 


1/2 


We  get  from  lemma  1.2 


1/2 


P(WM  <  s  '  )  =  M!fM(l)  , 


if  (fR)  is  the  y-Sheffer  sequence  for  D  with  roots  in  v,  where 


(3.1)  V1  =■  h”(i/M)  for  all  i  =  [c+(01)J  A  »1, . . . ,  |c+(02)_]  AM. 


and 


(3.2)  y±  -  h+(i/M)  for  all  i  «  [c_01)l+ . R"(02)l+-1  . 


1 


Ttie  following  short  tables  of  the  percentage  points  of  M1/2w„  are 

M 

computed  by  algoritym  A.l  and  by  the  outside  method  (1.5)  if  applicable. 
We  chose  always  0j_  -  0  -  l-02  for  6-0,  .01,  .05,  .1  and  .25.  Let 

P (*)  -  P(M1/2Wm  <  z)  . 

M 


We  consider  the  significance  probabilities  a  -  1-P(*a)  for 
a  *»  .1,  .05  and  .01.  Because  of  discontinuities,  these  levels  can 
not  always  be  attained.  If  the  absolute  difference  between  a  and 
l-p(za)  is  less  than  .000005  this  small  discontinuity  is  not  noted  in 
the  tables,  and  z^  is  rounded  to  4  digits  after  the  decimal  point.  If 


.000005  £  [a-l+P(za)|  <  .005  , 

and  a  is  greater (smaller)  than  1-P(z  ),  then  five  digits  are  given 

a 

and  a  bar  is  placed  under  (over)  the  last  digit.  This  last  digit  is  not 
rounded.  Decreasing  (increasing)  it  by  one  yields  a  probability  greater 
(smaller)  than  a.  Two  bars  indicate  an  absolute  difference  between 
.005  and  .013.  The  asymptotic  values  of  A. A.  Borokov  and  N.M.  Sycheva 
(1968,  Theorem  3A)  are  given  in  the  last  row  of  table  2-5. 

Table  1  is  a  confirmation  of  M.  Noe’s  (1972)  computations.  In  table 
2  and  3  the  results  of  P.L.  Canner's  (1975)  simulation  study  are  given  in 
parentheses.  In  table  4  and  5  the  rows  marked  by  F  contain  the  percen¬ 
tage  point 8  of  m*/2wm  as  the  tables  before.  The  rows  marked  by 
refer  to  the  correspondent  statistic  where  the  supremum  is  taken  over 
d/M  <  F  (x)  1-d/M.  The  integer  d  is  chosen  such  that  d/M  is 

closest  to  the  desired  0: 


(5-5) 


'1»J 

.[mo] 


if  M9- j_M0j  <  .5 


else  . 


The  F^-row  in  table  4  equals  for  M  ■  10  the  F-row  and  is  therefore 
omitted . 
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a=.l 


10  4.6146 


20  4.6423 


50  4.6631 


100  4.6719 


M 

a=-i 

a=.05 

a=  01 

10 

3.2900 

3-9829 

(3-33) 

6.03859 

(5.70) 

20 

3.3962 

4.04519 

(3.76) 

4.9094 

(4-77) 

5° 

3.2029 

3.55334 

(3.69) 

'4.5353 

(4.40) 

2-9218  3.4216  4.1705 

(3.33)  (4.00) 


2.9094  3.1831 
(3.07) 


2.8616 


1 100 j3. 06401 3.4379 


«M3.05  13-30 


100  2.8384  3.1417 


00  2.89  3.15 


M.  a=. 1  )  a=.05  a=. 01 


10. f  2.714813.13363.8203 


20 !f  2.713° -2.9830  3.72677' 

'■K  3-3938^4.0798  6.1725 

1  A  1 


50 ;F  2.7284j3.0071  3.6284 

Fx  2.9641 | 3. 3533  4.2777 


100 jF  2.7362i3.0120  3-5960 

! Fv  2.8562:3.1803  3.8839 


3.05  : 3  -  59 


Table  4:  8 


We  denote  by  W,,  ,,  the  two  sample  version  of  Ww: 

M,N  r  M 


W 

M,N 


sup 

0!<  Vx)^02 


|fx(x)-Fy(x) | 
[Fv(x)(l-Fv(x))]1/2 


(0^  =  a/(M+N)  and  0^  =  b/(JffN);  a  and  b  integer)  . 


Now  we  get  from  lemma  2.2 


N  „  „  1/2.  .  ,v.  ,,M+N. 

PWUM,NiS  >  '  £M<!!)/(  «  1  ' 


if  (f  )  is  the  y-3heffer  sequence  for  V  with  roots  in  V,  where 
n 


v.  =  rC^K)h'(i/H)1-i-l  and  ]i±  =  i_(tH-N)h+(i/M)J-i  , 


with  i  in  the  same  range  as  in  (3.1)  and  (3.2).  (See  (2.  5)  and  (2.6) 
for  v  and  y.)  The  following  tables  of  percentage  points  for  C^y^) ^ 
are  computed  using  only  algorithm  (2.^  ).  Discontinuities  occur  at  almost 
each  entry.  The  bars  are  set  following  the  same  rules  as  above,  but  only 
four  digits  are  given.  The  table  for  9=0  equals  the  table  for 
0  =  .01  and  is  therefore  omitted.  The  numbers  in  parentheses  are  taken 
from  P.L.  Canner's  (1975)  simulation  study  (computed  for  0=0).  For 
8  =  .05,  the  rows  M  =  N  =  10,  20  and  50  are  equal  to  those  in  table  6, 
and  are  omitted.  Instead,  we  demonstrate  the  effect  of  slightly  different, 
but  large  sample  sizes.  Again,  the  rows  are  marked  by  F^,  if  the  supremum 
is  taken  over  all  a'/M  £  Fx(x)  1  b’/M.  In  table  8  and  9  the  rows  are 
omitted  which  do  not  differ  from  table  6.  The  asymptotic  values  of  table 
2-5  may  be  used  for  comparison. 
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M=N  a=.l 


a=.05  |  a=.oi 


2.3W+1_  2.683 

(2.71)  (3-17) 
2.5819  2.7603  3-227^ 
(2.70)  (3-18) 


2.7OO7  2.9485  3.4299 
(2.93)  (3-44) 


2.791^  3-0249  3-5022 

(3.02)  (3-47) 


500  2.9441  3.I863  3-6694 


Table  6:  0  =  . 


a=.05  i  a=.0i 


Fy  2.7795  3-0082  j  3-5022 


Fv  2.7747 

3-01! 

Fy  2.7524 

2.18! 

Fv  2.7607 

3  011 

fy  2.7470 

2.99: 

Fy  2.7669 

3-02 

F.  2.7426 

i  2.99. 

500  500 


Table  7:  0 


M=W|  a=.i  '  a=.05  X  Q=.0l 


50  2.6667  .  2.8968  j  3-4299 


100  2.7003  i  2-9711  !  3-4720 


500  2.7415  3-0132  3-5468 


Table  8:  8 


M-N 

20 

2.3631 

50 

— 

2.4723 

dr=.05  a=.0l 


2-7357 

3  2965 

2-7580 

3.3131 

Table  9:  0  *»  .25 


For  applications  see  K.A.  Doksum  and  G.L.  Sievers  (1976):  "Plotting 
with  confidence:  Graphical  comparisons  of  two  populations." 
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.9515 

2.716326 

9  .9897  .9901 

3.259598 

9 

50 

50 

.8994 

.9013 

2.477696 

B 

.9490 

.9502 

2.735764 

9  .9895  .9901 

3.270847 

9 

50 

49 

.8953 

.9005 

2.480403 

8 

.9490 

.9502 

2.755113 

9  .9896  .9901 

3.268895 

9 

50 

48 

.8978 

.9002 

2.449985 

8 

.9496 

.9507 

2.746549 

9  .9897  .9904 

3.260541 

9 

50 

47 

.8971 

.9006 

2.467553 

8 

.9499 

.9510 

2.75722 7 

9  .9900  .9903 

3.304179 

9 

50 

46 

.8947 

.9007 

2.438465 

B 

.9495 

.9505 

2.744795 

9  .9894  .9901 

3.255123 

9 

50 

45 

.8989 

.9006 

2.491113 

B 

.9480 

.9510 

2.735209 

9  .9897  .9903 

3.264504 

9 

100 

100 

.8990 

.9002 

2.483005 

9 

.  9495 

,9502 

2.775122 

9  .9899  .9901 

3.313039 

9 

500 

500 

.9000 

.9001 

2.522724 

5 

.9500 

.9500 

2 . 802861 

9  .9900  .9900 

3.365530 

5 

The  purpose  of  this  appendix  is  to  summarize  the  algebraic  results 
which  we  used  in  the  preceding  chapters.  All  proofs  are  straightforward 
verifications  of  the  definitions.  Hence,  we  give  only  some  hints  and 
leave  the  details  to  the  reader.  A  more  general  approach  including 
Eulerian  polynomials  can  be  found  in  [ 11] .  The  "Finite  Operator  Calculus" 
of  G.-C.  Rota,  D.  Kahaner  and  A.  Odlyzko  [1^+  ]  is  the  fundament  of  the 
whole  theory. 

Let  P  be  the  algebra  of  polynomials  over  a  field.  K  with 
characteristic  zero.  In  our  rank  test  applications  K  always  equals 
7L  ,  for  the  order  tests  choose  K  =  1?  .  We  will  deal  with  linear 
operators  P  -♦  P  only,  and  omit  the  word  "linear"  in  the  sequel.  For 
all  a  e  K  the  shift  operator  is  denoted  by  Ea:  p(x)  p(x+a).  An 
operator  Q  on  ?  is  a  delta  operator,  if 

Q  is  shift  invariant:  QEa  =  EaQ  Va  e  K,  and 
Qx  is  a  non-zero  constant. 

The  derivative  operator  D  is  a  delta  operator  if  K  *  ]R  ,  and  the  following 
properties  show  how  Q  generalizes  D: 

(A.  1 )  Qa  =  0  for  every  constant  a  [14, p.687] 

(A.2)  deg(Qp)  =  deg(p)-l  for  each  p  e  P  with  deg(p)  >  1  [14, p.687)  . 

Hence,  the  kernel  of  Q  consists  only  of  the  constant  polynomials. 

A  sequence  of  polynomials  (sn)n  g  ^  is  a  Shef f er  sequence  for  Q,  if 


(A.3) 


sn  is  a  non-zero  constant 


(A*1*)  Qs^  =  for  all  n  >  1  . 

We  make  the  convention  8  “0  if  n  <  0.  For  instance,  (xn/n!)  is 

n 

a  Sheffer  sequence  for  D. 

Lemma  A.l;  If  (s^)  is  a  Sheffer  sequence  for  Q  then  deg(sn)  =  n. 
Proof;  (A.l)-(A.U)  ^ 

Lemma  A. 2:  If  ( s )  and  (tR)  are  both  Sheffer  sequences  for  Q 
with  the  property 

s  (v  )  =  t  (v  ) 
nv  n'  nv*n' 

for  a  given  sequence  (vn)  in  K,  then  the  two  sequences  are  equal. 
Proof;  Induction  over  n.  Use  ker(Q)  =  constant  functions  ^ 

(s  )  has  roots  in  v:  JL  -*  K ,  say,  if  s  (v  )  =  „ 

for  all  n  c  Eq  .  The  Sheffer  sequence  for  Q  with  roots  in  0  is 
called  the  basic  sequence  for  Q  and  always  denoted  by  (q_).  Obviously, 

(A.  5)  (xn/nl  )  is  the  basic  sequence  for  D  . 


It  is  easy  to  verify  that 


I 


(A. 6)  CX+n_1^nel}  *s  bas*c  se<luence  for  V  =  I-E"1  . 

More  examples  can  be  found  in  [1^] . 

Immediately  from  the  shift-invariance  follows:  If  (sn)  is  a 
Sheffer  sequence  for  Q  with  roots  in  v,  then  (E  sn)  is  a  Sheffer 
sequence  for  Q  with  roots  in  v-a. 

Deeper  than  all  the  other  results  in  this  appendix  is  the  following 

Lemma  A. 3:  If  v(n)  =  an+b  (a,b€K),  then 

sn(x)  :=  (x-an-b)(x-b)_1qn(x-b) 


defines  the  Sheffer  sequence  for  Q  with  roots  in  v. 

(For  n  =  0  we  have  to  define  —  ~  1.) 

Proof:  See  [1^,  p.  702]  ^ 

Now  we  come  to  a  representation  theorem  for  Sheffer  sequences  with 
roots  in 


v  (i ) 


r<P(i) 

.ci+d 


¥0  <  i  <  L 
Vi  >  L  , 


vhere  L  e  j  c,d£K  and  <p:  ]NQ  •*  R  arbitrary. 

Theorem  A. 1:  If  (s^)  is  the  Sheffer  sequence  for  Q  with  roots  in 
v  as  above,  then 


(A.7)  s  (x)  =  C  s,  (ci+d)(x-cn-d)(x-ci-d)“1qT,  ,  (x-ci-d)  Vn  e  H  . 

n  ito  1  n_i  0 
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Proof:  Check  recurrence  and  side  conditions,  using  lemma  A. 3 

_____  d 

Corollary  A.l:  (BinomiaiJCheor  em) .  If  (sn)  is  the  Sheffer-  and 
(q^)  the  basic  sequence  for  Q,  then 

n 

sn(x+y)  =  ^  si^^qLn-i^x^  Vn  e  ^0  * 

Proof:  Choose  c  =  0,  d  =  y  and  L  =  oo  in  (A.7)^ 

Avoiding  alternating  summation  in  (A. 7 ),  it  may  be  sometimes  preferable 
to  use  the  "outside  method"  (a  term,  introduced  by  J.L.  Hodges  (1957)): 

n  , 

(A.8)  sn(x)  =  rn(x)  -  Y.  (ci+d) (x-cn-d ) (x-ci -d )  q^Cx-ci-d) 

i=L+ 1 

where  (rn)  is  the  Sheffer  sequence  for  Q  with  roots  in  <p  (follows 

by  summation  over  all  i  =  0,...,n  in  (A. 7)). 

Repeated  use  of  (A. 7)  yields  a  representation  of  the  Sheffer 
sequence  (sn)  for  Q  with  roots  in  the  piecewise  affine  function 


-(i)  :=  iaj-Hjj  VLj  <  i  <  LJ+1  , 


where  -1  =  Ln  <  L,  <  •••  ,  each  L.  integer,  and  a.,b.  €  K  for  all 

U  X  J  u  J 


j  e  3Nq.  Uien  for  all  <  n  <  L^+1 


L1 

(A.9)  sn(x)  -  £  —  E  PJ(x)pJ_1(vJ(kJ))-...*p0(v1(k1))  , 


kj=0  k^ 
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if 


P,  (x)  = 


x-vi(ki^l) 
x-vi  (k± ) 


*k. 


1+1  i 


^(x-v.(k.))  , 


where  kQ  :=  0  and  :=  n-  Because  of  its  importance  we  explicitly 

write  down  the  special  case  of  (A. 9)  where 


ia+b 

• 

o 

tt 

5! 

ic+d 

Vi  >  L  . 

Then 


(a- io.'  v'x)  =  l  "ic+d^"  qi(ci+d~b)  Sfif  Vi(x"ic-d) 


i=(. 


For  n  <  L,  the  r.h.s.  equals  (x-an-b ) (x-b )-1qn(x-b )  by  lemma  A-3 


Now  we  assume  that  K  is  completely  ordered. 

be  a  non  decreasing  function  and  (t  . )  .  _T 

x  n,i'n,i  e  HQ 

in  P  with  the  properties 


Let  p:  B0+K 
be  a  double  sequence 


(A.  11) 


Vi^j.)  =  *  ^(hj)  VO  <  i  <  r(n)  :=  min{ra  e  3NQ| p(m)=p(n)), 

=0  Vi  >  r(n)  . 


Define  an  associated  sequence 


(fn>  to  by 


(A.  12 )  fn(x)  tn^(x)  V4i-1  <  x  <  ^  (p^  :=  -«)  . 

We  call  (fn)  a  p-Sheffer  sequence,  if  (t^ r(n0)n  e  *  Is  •  Sheffet 
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1 


sequence  for  all  m  o  Kg  .  From  corollary  A.l  we  get  a  first  represen¬ 


tation  of  f  (x) : 

n 


(A 


n 

.13)  V*)  =kC  fk(y\.k(x-y)  if  x,y€[Mi_r4jL]  . 


If  (fn)  has  roots  in  v,  i.e. 


(A. 14) 


W  '  50,n  V"  '  ®o  ’ 


then  any  value  fR(z)  can  be  computed  from  (A. 13)  by  stepping  through  all 
the  intervals  until  z  is  enclosed.  We  give  only  a  brief 

description  of  this  trivial  algorithm: 


Algorithm  A.l.  Assume  ^r(j)  (V*j)  >  •  •  •  >f  )  are  already  computed  such 
j  £  n  at 
a)  If  v. 


that  j  £  n  and  ^i+1  > 


'i+1  <  Uj+1  then  define  x  :=  v±+1,  y  :=  p 
and  compute  fr^  j  (v±+1>  , . . .  .f^v^)  from  (A. 13).  Of  course, 

fi+l(Vi+l)  =  Therefore,  the  i-index  increased  by  one,  and  it  increases 
again  if  Hes  also  in  the  same  interval  (define  x  ■>  and 


y  -  vi+1)-  Finally  a  k  is  reached  such  that  y(j)  <  <  y^+1  <  vfc+1 

(the  case  =  ^j+1  is  to  c^e  reader).  Then  choose 

x  :=  y.,,  ,  y  :=  V,  ,  and  compute  f  ...  (u  f  \ 

3+1  k  *  r(j)  ^j+1 . V’j+r 

from  (A. 13).  Now  we  are  in  the  same  situation  as  in  the  beginning. 

b)  If  v±+1  >  yj+1  >  then  define  x  :=  yj+1,  y  Pj 

compute  ^r(j)  *  *  *  ’^k^j+l^  from  (A. 13).  Again,  we  are 

in  the  same  situation  as  in  the  beginning. 

c)  If  y^  ■  y^+1  increase  j  by  one. 

In  special  cases  this  algorithm  can  be  simplified. 


and 


53 


i  _ 


! 
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A  one  dimensional  recursion  can  be  obtained  from 

Theorem  A.  2:  Let  (fR)  be  a  y-Sheffer  sequence  for  Q  (with  basic 
sequence  (q^).  If  (fn)  is  associated  to  (tR  i>  then 

n 

(A. 15)  tR  i(x)  =  l  fk^k^qn-k^x_yk^  for  a11  n  E  W0  and  i“°»*--»n  • 

k=  i 

Proof :  Verify  side  conditions  (A. 11). 

See  (26,  theorem  4.1]  for  a  general  version  of  this  theorem.  The 
announced  one  dimensional  recursion  follows,  when  we  write  (A.  15)  as 

(A- 16)  fn^x)  =  fk^k^qn-k^X'Wk^  for  a11  n  e  , 

where  the  summation  runs  over  all  k  such  that  >  x.  Thus,  a  system 
of  equations  for  the  unknown  fj^H^)  is  obtained,  if  only  one  value 
f  (v  )  with  V  <  y  is  known  for  each  n.  By  Cramer's  rule,  f  (u  ) 
can  be  expressed  as  a  determinant: 

if  (fjj)  is  a  y-Sheffer-  and  (qn)  the  basic  sequence 

for  Q,  then 


fnV 

,n+l  »  where 

CT  ,  = 

f  Vjtt'i-l-Mj.i).) 

VJ  « 

if  i 

1  fi-l  ^Vi-l) 

if  i  =  n+1  t 

xor  any  v  <  y.  If,  iQ  addition,  (fj  hM  ^  ^  ^ 
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